Derivative moments in turbulent shear flows 
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We propose a generalized perspective on the behavior of high-order derivative moments in tur- 
bulent shear flows by taking account of the roles of small-scale intermittency and mean shear, in 
addition to the Reynolds number. Two asymptotic regimes are discussed with respect to shear ef- 
fects. By these means, some existing disagreements on the Reynolds number dependence of derivative 
moments can be explained. That odd-order moments of transverse velocity derivatives tend not van- 
ish as expected from elementary scaling considerations does not necessarily imply that small-scale 
anisotropy persists at all Reynolds numbers. 
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I. INTRODUCTION 

The postulate of local isotropy 0] implies an invariance 
with respect to spatial rotations of the statistical proper- 
ties of small scales of turbulence. Even though the large 
scales are anisotropic in all practical flows, it is thought 
that the small scales at high Reynolds numbers are 
shielded from anisotropy because of their separation by a 
wide range of intermediate scales. At any finite Reynolds 
number, some residual effects of small-scale anisotropy 
may linger, but all proper measures of anisotropy are ex- 
pected to decrease rapidly with Reynolds number. An 
understanding of the rate at which small scales tend to- 
wards isotropy is a basic building block in turbulence 
theory. 

A particularly appealing manner of generating large- 
scale anisotropy is by a homogeneous shear characterized 
by a constant shear rate S = dU /dy, where U (y) is the 
mean velocity in the streamwise direction x, and y is the 
direction of the shear. During the last few years, nearly 
homogeneous shear flows, both experimental |^,[|,fl| and 
numerical, [^fl^J^J have examined the rate at which lo- 
cal isotropy is recovered with respect to the Taylor mi- 
croscale Reynolds number, R\. The notation is standard: 
R\ = u'X/v, u is the velocity fluctuation in the longitu- 
dinal direction x, vl = \J (u 2 ), A 2 = (u 2 )/((du/dx) 2 }, v 
is the kinematic viscosity and (•) denotes a suitably de- 
fined average. The discussion has often been focused on 
the behavior of normalized odd moments of transverse 
velocity derivatives defined as 

((du/dy) 2n+1 ) 
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where n is a positive integer. The velocity derivatives 
are small-scale quantities, and symmetry considerations 
of local isotropy demand that the odd moments of trans- 
verse velocity derivatives be zero. In practice, they should 



decrease with R\ relatively rapidly. Though the postu- 
late of local isotropy does not by itself stipulate this rate, 
simple estimates can be made by retaining the spirit of 
local isotropy and making further assumptions. Let us 
assume that the non-dimensional moments Mm+\ de- 
pend on some power p of the shear, which is the ulti- 
mate source of anisotropy in homogeneous shear flows, 
and non-dimensionalize the shear dependence by a time 
scale formed by the energy dissipation rate e and fluid vis- 
cosity v; this non-dimensionalization accords with Kol- 
mogorov's first hypothesis that small-scale properties de- 
pend solely on e and v. We may then write 
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(du/dy) = S p f(e,v) 
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where it is further assumed that S = 0(u'/L) and L is 
the large scale of turbulence. Lumley Q considered a lin- 
ear dependence on the shear (i.e., p = 1). This yields an 
inverse power of R\ for the decay of all odd moments. 
The choice p = 1 accounts for the dependence of the sign 
of the odd moments on the sign of S in a simple manner. 

The existing experimental and numerical data on the 
skewness of the transverse velocity derivative are col- 
lected in Fig. 1. Data from any given source seem to fit 
a power-law of the form 



M 3 (du/dy) ~ R: 



(3) 



with to differing from one set of data to another. Though 
the power-law is not a particularly good fit for the totality 
of the data, the average roll-off seems to be less steep than 
the inverse power discussed above. Atmospheric data at 
much higher Reynolds numbers |} 10 are consistent with 
this slower rate of decay. 

The situation with respect to the hyper ksewness, M§, 
is as follows. Two independent laboratory experiments in 
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nearly homogeneous shear flows |3j||] draw different con- 
clusions on the R\ behavior of M5. On the one hand, 
Shen and Warhaft find no dependence on R\ in the 
range between fO 2 and fO 3 . On the other hand, Ferchichi 
and Tavoularis || regard their data to be essentially con- 
sistent with expectations of local isotropy. (For one per- 
spective on this difference, see Warhaft and Shen [[!)). A 
collection of all known data is given in Fig. 2. The overall 
impression from the figure is that, while there is probably 
a weak decreasing trend for R\ > 300, the hyperskewness 
does not diminish perceptibly even when R\ is as high 
as f 000. Shen and Warhaft Q measured the seventh nor- 
malized moment of the velocity derivative and found that 
it increased with R\ instead of decreasing. 

Such findings have been interpreted (e.g., Refs. 4 and 
5) to mean that, in the presence of large-scale shear, 
small-scale anisotropy persists to very high Reynolds 
numbers. The intent of this paper is to clarify, at least 
partially, the simultaneous role played by shear, intermit- 
tency and the Reynolds number — all of which have an 
impact on trends displayed by odd derivative moments. 

We now consider in Sec. II the issue of intermittency 
versus anisotropy. In Sec. Ill, we highlight the effects of 
shear by discussing the limiting cases of large shear and 
local isotropy, and argue that R\ determines the state of 
the flow only partially. In Sec. IV the derivative skew- 
ness data are plotted in the plane spanned by the non- 
dimensional shear parameter S* and R\, where we use 
the definition 

S* = Su' 2 /e. (4) 

This broader perspective may help resolve some seem- 
ingly contradictory claims on the recovery of isotropy of 
small scales. 



II. INTERMITTENCY VERSUS ANISOTROPY 

A large number of measurements have shown convinc- 
ingly that high-order even moments of small-scale fea- 
tures of turbulence increase with R\. Consider the lon- 
gitudinal velocity derivative du/dx. The product of its 
second moment and fluid viscosity is essentially the en- 
ergy dissipation, which is known to remain independent 
of R\ when R\ exceeds some moderately high value (see 
Refs. 14,15). On the other hand, all high-order moments 
of longitudinal velocity gradients grow with R\ (see, for 
example, the compilation in Ref. 16 of the data on the 
flatness factor of du/dx). Sixth and higher order mo- 
ments increase at increasingly faster rates with R\ . These 
growths are attributed to the intermittency of small- 
scale turbulence. At the present level of our understand- 
ing, intermittency is independent of anisotropy effects. 
Therefore, just as the growth of high-order even moments 



with increasing R\ is unrelated to anisotropy, it is legiti- 
mate to ask if, at least in part, the slower-than-expected 
decay — or even modest growth — of odd moments, may 
be related to intermittency. 

To separate intermittency effects from those of 
anisotropy, at least in some approximate way, it is useful 
to plot the ratio —M2 n +i{du/dy)/M2 n +i(du/dx). It is 
plausible to assume that the i?A-growth due to intermit- 
tency effects is the same for the moments of du/dx and 
for the moments of du/dy, the intermittency effects of 
odd moments of du/dy are cancelled in these ratios by 
those of du/dx. Though this not a rigorous statement, it 
is useful to see the outcome. Figure 3 shows the results. 
It is clear, despite large scatter, that all the moments 
show a tendency to diminish with Reynolds number. The 
odd- moments in Fig. 3 are normalized by powers of their 
respective variances. It would have been desirable to plot 
ratios of unnormalized, or "bare" , moments of the two 
derivatives, but Ref. 4 does not include those data. In 
any case, this should not make much difference because 
((du/dx) 2 } /((du/dy) 2 } is essentially a constant at high 
Reynolds numbers. 

This same issue can be rephrased and reexamined in 
a somewhat different light. When we consider the mo- 
ments such as skewness and hyperskewness, we usually 
normalize them by the appropriate power of the variance 
of the variable. This is perfectly reasonable for Gaussian 
or near-Gaussian variables, but not so for intermittent 
quantities with highly stretched tails. Perhaps a more 
reasonable alternative is to consider how an odd moment 
of a certain order varies with respect to the even moment 
just below or just above, or the geometric mean of those 
just below and just above. We illustrate the results of 
this consideration for the third, fifth and seventh order 
moments of du/dy (= x) in Fig. 4. The lack of data on 
the eight moment of du / dy makes the analysis of the sev- 
enth moment incomplete. Nevertheless, it is clear that all 
these alternative ways of normalization show substantial 
decay. It is hard to be precise about the rates of decay, 
partly because of the large scatter and partly because 
the incomplete manner in which the seventh moments 
have been analyzed, but it is conceivable that increasingly 
high-order moments, within a given normalizing scheme, 
decay more slowly. At the least, a careful discussion of the 
restoration of anisotropy requires the proper inclusion of 
intermittency effects. This is our first point. 



III. SHEAR EFFECTS IN TWO LIMITING CASES 

It is reasonable to suppose that, to a first approxima- 
tion, the mean shear S and the viscosity v determine the 
gross state of the flow. Expressing time units in terms of 
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S 1 , length units in terms of the integral length scale L, 
and the mean profile as U = Sye x , we get 

— + (v • V)v + Sy— + Sve x = -Vp + ^V 2 v + f , 

V ■ v = , 

where v = (u, v, w) and e x is unit vector in the stream- 
wise direction x. The volume forcing is denoted by f. 
The two parameters may be expected to set the steady 
state fluctuation level and energy dissipation rate. It also 
follows that the derived parameters R\ and S* adjust 
themselves dynamically, in ways that are understood only 
partially, to the imposed values of v and S . We then ex- 
pect 



least, we have to possess some knowledge of the other pa- 
rameters influencing the state of the flow before drawing 
firm conclusions on the recovery of isotropy. 

To keep matters simple, we will focus below on homo- 
geneously sheared flows. Because the initial conditions 
(6) are not known in all quantifiable details, we shall tenta- 
tively stipulate a simple generalization of Eq. (1) in the 
form 



R\ = g\{v,S) = gi(Re s ) 
S* = 92{v,S) = 92{Re a ) 



(7) 



where g\ and gi are unknown functions of their argu- 
ments. 

All the homogeneous shear experiments to be discussed 
below have been done in air. This fixes the viscosity to be 
approximately constant, so we can test the dependence 
of S* and R\ on the shear rate S. The relevant plots 
(Figs. 5 (a) and (b)) show no obvious trend but only 
large scatter. This scatter may be related in part to the 
fact that P/e ^ 1 in some of the experiments, leading 
to nonstationarity. Here P is the production of turbulent 
kinetic energy, defined as P = — (uv)S for homogeneous 
shear. In part, it demonstrates that the flow might de- 
pend additionally on initial conditions, (v 2 ) t= 0j or the 
type of driving of small-scale fluctuations characterized 
globally by an input energy. From Eq. (JsJ) , the latter is 
given by 



We then get more complex relations such as 

R\ = gi{v, S, t in ) = <7i(-Re s ,e in ) , 
S* = g%(y, S,t in ) =gi{Re s ,t in ). 



(8) 



(9) 



If we add to Figs. 5 the atmospheric data from Refs. 9 
and 10, the situation becomes even more complex. How- 
ever, it is likely that, in such inhomogeneous flows, one 
has to take into account secondary factors such as convec- 
tive effects (though conditions in which they are modest 
can always be chosen carefully). In laboratory experi- 
ments, secondary effects might arise from the use of spe- 
cific passive or active grids for the generation of turbu- 
lence. Further differences can arise when measuring at 
a fixed point instead of following the downstream evolu- 
tion. This discussion merely underlines the inadequacy of 
R\ as the sole indicator of the state of the flow. At the 



S 3 = f(R x ,S*), 



(10) 



and regard other effects as superimposed "noise" p5[ . 
If so, we should investigate the behavior of the deriva- 
tive moments by keeping one of the two parameters fixed 
while varying the other, for example by fixing S* and 
varying R\. This is the topic of the next section. How- 
ever, it is useful to preface this consideration by exam- 
ining two limiting behaviors in which some inequalities 
between R\ and S* can be established. 



A. Large shear case 

Consider the case of large shear rate for which the cou- 
pling of the mean shear to the small-scale flow dominates. 
In the rapid distortion limit, namely S —> oo, Eq. (|5|) 
becomes linear because the viscous term as well as the 
(v • V)v term can be dropped, so that any shear rate de- 
pendence can be eliminated by the rescaling of the vari- 
ables, e.g. t — > St. 

Our dimensional estimates are related to large but fi- 
nite shear rates. For this case, the term representing the 
coupling of the turbulent velocity component to the mean 
shear is important and large in comparison to the non- 
linear advective transport. Our situation corresponds to 
the case in which 
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dvj_ 

' dxi 



(11) 



For the homogeneous shear flow, we get for the left-hand- 
side of this equation 



dU 
dy 



= \v\S- 



(v 2 ) 1 / 2 S<(u 2 ) 1 / 2 S. 



(12) 



In the last step above, we have used the fact that the 
root-mean-square velocity in the streamwise direction is 
larger than that in the shear direction, as has been found 
in all numerical and physical experiments. The term on 
the right hand side can be estimated roughly as (u 2 )/£ 
where the scale £ is characteristic of turbulent velocity 
gradients, and can therefore be assumed to be of the or- 
der of the Taylor microscale, A. We then require 



(13) 
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With e = c e u ri /L and L/X = Re^ 2 /\/l0, and the con- 
stant c e ~ 0(1) for sufficiently large R\, we get 



S* y/3 
Rx V200 



(14) 



In reality, c e depends weakly on S"*. For example, Sreeni- 
vasan [£3| has examined the data and concluded that 
c e » c G exp(— aS 1 *), where a ~ 0.03, is a good empirical 
fit. Since this dependence is quite weak, we have taken 
c e to be a constant for simplicity. Presumably, if Eq. (14) 
holds, the effects of shear will always be felt no matter 
how high the Reynolds number. 



B. Local isotropy limit 

At the other extreme is the case in which local isotropy 
can be expected, a priori, to prevail. A suitable criterion 
(see Corrsin J24J) is that a sufficiently large separation 
should exist between the shear time scale, S , and the 
Kolmogorov time scale, t v — (is/e) 1 / 2 . This can be writ- 
ten as 



or as 



St„ « 1 , 



ST v = S*JFAr-Re- L 1/2 . 



(15) 



(16) 



With c e ~ O(l) and R x = (20i?e L /3) 1 / 2 we obtain 

0.35. (17) 



S* 

77" < 17= 
Rx \/200 



The implication is that local isotropy prevails for all 
S*/R\ substantially smaller than 0.35. For all other con- 
ditions, one should expect that the magnitude of S* will 
play some role in determining how high an Rx is required 
for local isotropy to prevail. This explicit dependence on 
shear has been noted for passive scalars by Sreenivasan 
and Tavoularis; pq| see their figures 2 and 3. 



the turbulent energy and on the energy dissipation esti- 
mate from the local isotropy formula, e = 15v((du/dx) 2 ). 
The error made in this estimate for the energy dissipa- 
tion depends [|2j| on the magnitudes of shear and R\, 
but it appears to be a reasonable approximation for the 
present conditions. We have recalculated for all numeri- 
cal data the energy dissipation rate as in experiments. It 
is clear from the figure that there is no simple correlation 
between the two parameters S* and R\. 

We have replotted the same data in Fig. 7. Different 
symbols correspond now to different magnitudes of the 
skewness, as indicated in the legend, and not to different 
experiments. Superimposed are islands in shades of grey 
to obtain a rough idea of the surface plot of / (R \ , S* ) . We 
used an interpolation routine with local thin plate splines 
that can reconstruct a surface from scattered data. |2S| ] 
The data are sparse in most of the regions of the plane 
so any surface-fitting routine will introduce some pecu- 
liarities. For reasons explained above, we can expect that 
the derivative skewness will become small in the local 
isotropy limit (lower right corner) and that the values 
will grow above 0.8 in the large shear case (upper left 
corner). Although the latter fact is not reflected by the 
surface fit because the data points are absent there and 
almost all data points are in the intermediate region, we 
think that the surface plot is not unreasonable. 

This plot offers additional perspectives. For instance, 
Fig. 1 is just a projection of the data on to the R\ axis 
and masks the fundamental effect of the applied shear 
(among other effects). We have marked in this figure 
the trends of the data sets of Warhaft and co-workers 
(parabolic solid line labelled by GW/SW) and of Fer- 
chichi and Tavoularis (straight line labelled by FT), re- 
spectively. They show that the two experiments follow 
different paths: while the Ferchichi- Tavoularis data run 
directly down the 'mountain range', the Shen- Warhaft 
data seem to run in a kind of spiral around the 'moun- 
tain', presumably resulting in a slower decay of the skew- 
ness when projected on to the R\ axis. Finally, the two 
limiting regimes of Sec. Ill are also plotted. The local 
isotropy limit is not reached for any of the data. The 
large shear case is reached for very large S* . A threshold 
S* grows when R\ increases; alternatively, local isotropy 
requires larger Rx if the shear parameter is larger. 



IV. THE R x -S* PHASE DIAGRAM 

We now plot in Fig. 6 all available data on the skew- 
ness of the transverse derivative, du/dy, on a phase plane 
consisting of S* and Rx- The conventional normaliza- 
tion factors in the definition of S* are the total turbulent 
energy and its dissipation rate. This can be done quite 
readily for the numerical data, but experiments usually 
provide information only on the streamwise component of 



V. CONCLUDING REMARKS 

Considerable attention has been paid recently to the 
fact that in shear flows the moments of transverse ve- 
locity derivatives do not vanish with Reynolds number 
as fast as was expected. We have introduced two con- 
siderations for interpreting these observations, invoking 
small-scale intermittency and the magnitude of the shear 
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parameter. These two effects work in combination with 
the Reynolds numbers in determining the magnitudes of 
odd moments of velocity derivatives. The fact that the 
odd moments, when normalized by an appropriate power 
of the variance (a procedure steeped in studies of Gaus- 
sian processes) , decay more slowly than expected should 
not be considered a priori as incontrovertible evidence 
against local isotropy. We believe that the broader per- 
spective of this paper may explain some seeming contra- 
dictions that exist in the literature. 

The conclusions we draw in this paper would be more 
definitive if the data spanned much higher Reynolds num- 
ber range. This can be done with adequate resolution 
only in atmospheric flows at present. The existing mea- 
surements are in conformity with the discussion here, but 
it is difficult to be definitive because of the usual prob- 
lems that often exist in atmospheric measurements. On 
the other hand, it seems quite feasible in numerical simu- 
lations to fix the Reynolds number and vary the shear pa- 
rameter, even though the Reynolds number range may be 
limited. Such a study will tell us more about the restora- 
tion (or otherwise) of small-scale isotropy. 

In the recent past, various efforts have been made to 
understand the effects of anisotropy through an SO (3) 
decomposition of structure functions and other tensorial 
objects (see Ref. 29 for the basic theoretical idea and 
Rcfs. 30-32 for implementations of the idea and further 
references). The method offers a transparent way of de- 
termining the degree of anisotropy in turbulent statistics. 
The relation between that effort and the present global 
picture needs to be explored. 
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FIG. 1. The dependence of the skewness of the transverse velocity derivative on the Taylor microscale Reynolds number. The 
solid line indicates the R\ scaling expected to hold on the basis of Lumley's dimensional considerations (see text). Notation in 
the legend stands for: PU, Pumir [5]; YE, Yeung, for details of whose shear flow simulations, see [12]; RM, Rogers et al. [13]; 
SE, Schumacher and Eckhardt [6,7]; GW, Garg and Warhaft [2]; SW, Shen and Warhaft [4]; FT, Ferchichi and Tavoularis [3]. 
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FIG. 2. The dependence of the hyperskewness of the transverse velocity derivative on the Taylor microscale Reynolds number. 
The notation in the legend is the same as in Fig. 1. 
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FIG. 3. The ratio of the normalized moments M2 n +i(du/dy) and M2 n +i(du/dx) for n = 1, 2 and 3 versus the Taylor 
microscale Reynolds number. Quantities M 2n +i are denned in Eq. (1). The notation in the legend is the same as in Fig. 1. Here 
and elsewhere, data from [4] have been read from published graphs. 
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FIG. 4. Third, fifth, and seventh order transverse derivative moments with different normalizations, plotted against R\. 
Squares are normalized by the moment of preceding even order, circles by the succeeding even order. Triangles are normalized 
by the geometric mean of both orders. Variable x in the legend stands for du/dy. The data are from Yeung [12], Schu- 
macher/Eckhardt [6,7] and Shen/Warhaft [4]. Solid lines with attached numbers indicate fits to possible algebraic power laws. 
Each symbol corresponds to a different normalization, and includes all the data just cited. 
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FIG. 5. (a) The dependence of the shear parameter S* on the shear rate S for fixed kinematic viscosity, v — v a i r . (b) 
Dependence of the Taylor-microscale Reynolds number _Ra on S for fixed v. Unfilled symbols are the same as in Fig. 1. 
Additional shear flow experiments have been included, though they did not focus on derivative moments explicitly (filled 
symbols). The additions are: RO, Rose [17]; ML, Mulhearn and Luxton [18]; HC, Harris, Graham, and Corrsin [19]; TC, 
Tavoularis and Corrsin [20]; TK, Tavoularis and Karnik [21]; SV, Saddoughi and Veeravalli [22]. 
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FIG. 6. Operating points of the homogeneous shear flows in the R\-S* plane. All symbols are the same as in Figs. 1 and 5. 
Data points LK stand for Lee et al. [26], corresponding to additional numerical simulations at very high shear rates. 
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FIG. 7. Derivative skewness and its dependence of R\ and S*. The solid straight line is for the large shear limit S*/R\ = 1, 
the dashed lines are for local isotropy limit S* / R\ = 0.003 and 0.001 (left to right). Underlying grey scales result from a surface 
fit. Trends of the data from Ferchichi and Tavoularis (FT) and Warhaft and co-workers (GW/SW) are indicated by solid lines. 
Different symbols indicate different ranges of the skewness, and may represent data from the same source. The data point below 
S* = 3 is taken from the numerical experiments of Rogers et al. [13] 
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